Energy transfer dynamics and thermalization of two oscillators interacting via chaos 
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We consider the classical dynamics of two particles moving in harmonic potential wells and interacting with 
the same external environment He, consisting of A' non-interacting chaotic systems. The parameters are set so 
that when either particle is separately placed in contact with the environment, a dissipative behavior is observed. 
When both particles are simultaneously in contact with He an indirect coupling between them is observed only 
if the particles are in near resonance. We study the equilibrium properties of the system considering ensemble 
averages for the case N = I and single trajectory dynamics for A' large. In both cases, the particles and the 
environment reach an equilibrium configuration at long times, but only for large A' a temperature can be assigned 
to the system. 



INTRODUCTION 

Understanding dissipation from a microscopic point of 
view has become important to several areas of physics, es- 
pecially if quantum phenomena are relevant. In these cases 
a globally conservative approach for the system plus its en- 
vironment is highly desirable, allowing a direct quantum me- 
chanical description. 

The simplest, and perhaps most natural, way to model the 
environment is to use an infinite set of harmonic oscillators, 
representing the normal modes of a general system in equilib- 
rium weakly perturbed by the system of interest Iilii2j. The 
spectral function, which is related to the distribution of fre- 
quencies of the normal modes, can be chosen to model sev- 
eral types of thermal baths OKS). Other representations of the 
environment have also been explored, from spin systems ||6l 
(or two-level atoms) to chaotic systems [7-11 1. The latter is 
particularly important to model coupling to small external sys- 
tems where the chaotic nature of the trajectories compensates 
for the small number of degrees of freedom in the decay of 
correlation functions. However, chaos alone does not suffice 
to simulate a thermal bath, since small numbers of degrees 
of freedom always leave a strong signature in the dynamics 
through large fluctuations in the observables of the system of 
interest lfT2l . These fluctuations can be washed out by aver- 
aging over several realizations of the dynamics lfT3l . 

More recently, the interest has shifted from a single parti- 
cle interacting with the environment to two particles indepen- 
dently connected to it lll4HT6l . allowing the study of interac- 
tions mediated by the environment. One important case is that 
of two entangled particles subjected to dissipation and deco- 
herence. In many cases, the total Hamiltonian is symmetric 
by the exchange of the particles, although they may still be 
considered distinguishable in some applications [14J. 

In this paper we consider the classical dynamics of two 
particles moving in harmonic potentials linearly coupled to 
a finite chaotic environment. We study the system behavior as 
a function of the frequency of the oscillators. We find that the 
behavior of the system changes radically when the oscillators 
are in resonance, which is a necessary condition if the particles 
are identical, and that even small deviations from this sym- 
metric state changes dramatically the equilibrium properties 



of the system. In particular, we show that, when in resonance, 
the particles may exchange energy through the environment 
while their energies dissipate. Moreover, the energy stored in 
their relative motion is conserved. 

We model the environment as a set of independent quar- 
tic systems (QS), each with 2 degrees of freedom. The QS 
has a single parameter that controls the degree of chaos in the 
dynamics. The particles are represented by two harmonic os- 
cillators independently coupled to the set of QS's. As study 
cases we consider the systems with N — I and A^ = 100, for 
which we study the statistical properties of the energy dis- 
tributions of the environment and of the oscillators. In both 
cases, these distributions reach an asymptotic equilibrium, but 
attempts to define a temperature for the system using two ba- 
sic definitions of entropy, applicable to small systems, suc- 
ceed only in the case A^ 100. The results of the simulations 
are then interpreted in the light of the linear response theory. 



MODEL 

We consider two harmonic oscillators interacting with an 
environment composed by a collection of independent quartic 
systems. The Hamiltonian is 



H = Hi+H2+He + XnHi, 



where Xn = X / VN and 
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The total energy is conserved and the two oscillators interact 
only via the environment. 

In our simulations, we used A = 0.01. The parameter a 
in Hqs controls the dynamical regime of the quartic systems 
in the environment, ranging from integrable (for the special 
values fl = 1.0 and a = 0.33 ifTTl ) to chaotic (a — > 0). In 
the work, we used a = 0.01 or a = 0.1, which correspond 
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to regimes where the QS is mostly chaotic. For a = 0.1 the 
largest Lyapunov exponent of an isolated QS is Xi = 0.166 
for energy Eqs = 0.1. For a = 0.01 we obtained Xl = 0.121 
for Eqs = 0.01. Because of the scaling properties of the QS, 
the Lyapunov exponent at other energies can be calculated us- 

ingAL(£) = (£/£o)'/'Ai(£o). 

As the number of degrees of freedom is finite, the parameter 
A' assumes a prominent role and we investigate its influence in 
the dynamic behavior of the oscillators. It has been recently 
shown [12 1 that, for sufficiently large, such finite chaotic 
environment can simulate the action of an infinite thermal 
reservoir. Here the environment also acts as a medium con- 
necting the two harmonic oscillators. 

In order to highlight the interaction between the oscillators, 
one of them is initialized with energy Ei > while the other 
is set at rest with E2 = 0. For the environment we define the 
initial conditions using the pseudo canonical distribution l,12J 
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where the energy Eq^ of each QS is randomly chosen from the 
exponential probability distribution exp {—E /Eqs) /Eqs- The 
value of Eqs plays the role of an initial "temperature" for the 
environment. In the special case of = 1 the energy Eqs is 
fixed to Eqs- 

In what follows we will study the system ([T]) for only two 
relevant environment sizes: the "microscopic" (A^ = 1) and 
"macroscopic" (A^ — 100) cases. As pointed out in lfT2ll . the 
dynamical behavior of the system becomes A^-independent for 
A^ sufficiently large. For the parameters values used in this 
paper, the large A^ limit is already reached for A^ = 100. The 
case A^ 1 is a natural extension of the work presented in 
ifTSl . where the authors considered the interaction between a 
harmonic oscillator and a single quartic system using ensem- 
ble averages. Thence we choose the same set of parameters 
as in lfT3l in our simulations, allowing us to verify the im- 
plications of adding a second oscillator to the system. For 
A^ = 100, on the other hand, we compare our results with the 
work presented in IIT2I . which treated the dynamics of a single 
oscillator interacting with large chaotic environments. 

In the microscopic case, observables related to the har- 
monic oscillators, like the energies Hi or H2, exhibit large 
fluctuations when coupled to Hqs- These fluctuations can 
only be washed out by averaging over large ensembles of re- 
alizations of the dynamics. In the macroscopic case this is not 
necessary and the results obtained from a single realization 
of the dynamics are already representative of the average be- 
havior In this case we can speak of an "effective dynamics", 
where no averages are performed. 

MEASURES OF TEMPERATURE 

Temperature is a central property in the description of equi- 
librium and is properly defined only in the thermodynamic 



limit of very large systems. Since the environment defined in 
Eq. (|3]l is far from this limit for N — I and A^ — 100, different 
possibilities arise. One natural definition comes from the 
equipartition theorem IfTSlI 
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where z„ denotes the coordinates or momenta of the system 
and 5i„n is the Kronecker delta. The constant Te is identified 
with ksT when the system is in contact with a thermal reser- 
voir The subscript in Te emphasizes the explicit use of the 
equipartition theorem. Eqs. ([T]i and (|6| predict that the frac- 
tion of the total energy within each subsystem in equilibrium 
should be (Ei) = {E2) — Te = '^■/^{Eqs)- We can also write 
the system's temperature as a function of the number of quar- 
tic systems as 

Et = £i(0) +£2(0) +4s(0) =2te + ^Nte, (7) 

in which Eg^lO) is the environment's total energy at f = 0. 

Eq. Q can also be derived from the thermodynamic rela- 
tion |fT9ll 



dlnF 
dEr 



where 



r(£r) = J ®{ET-H)dpdq 
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and the integral is taken over the entire phase space of the 
system. For the Hamiltonian in Eq. ([T]) we obtain 



T^cE- 



(10) 



where Nf ~ 2N is the number of degrees of freedom of the 
environment and c represents a constant depending on the sys- 
tem parameters. 
Notice that 



InF = Inc + (3A^//4 + 2) InEr 



(11) 



plays the role of entropy. The usual entropy, on the other hand, 
is given by 



S{E) = In 



5{E — H)dpdq 



= ln 
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In c + (3A^/-/4 + 1 ) ln£ + In (3A^/-/4 + 2) . (12) 



Defining ' = dS{ET) / dEj we obtain 



Et = Tj+ -A^T5, 



(13) 



which agrees with (|7| in the limit of large A^. 

The two temperatures, Te and T5, can be calculated numer- 
ically in a number of ways. According to |20l, Ts may be 
calculated from the dynamics using the expression L20i,21J 
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Figure 1. (Color online) Short- and long-time energy dynamics of the oscillators at various values of k = (black), 1 (purple), 6 (green), 11 
(orange), 16 (cyan) for energies of Eqs{0) =0.1 (panels (a) and (b)) and Eqs{0) = 5 (panels (c) and (d)). The averages were computed using 
an ensemble of 40000 initial conditions. Colors change gradually from light gray to black as the k value decreases, if printed in black and 
white. 



where t^'/^b the so-called Rugh's temperature and 

with V denoting the gradient in phase space. This formula as- 
sumes that the energy, H = E, is the only isolating integral. 
On the other hand, Te can be estimated as a dynamical aver- 
age, as in the left hand side of Eq. (|6]l. We use the superscript 
d to emphasize that these temperatures are obtained from the 
dynamics, i.e., from the trajectories. 

Finally, the temperature can be calculated from fitting the 
energy distribution p{E), which is the probability of finding 
one QS with energy E when the system is in equilibrium, 
as the Boltzmann exponential pb{E) ^ exp(— E/tb). An- 
other possibility is to numerically calculate the distribution 
of momentum and fit it with the Maxwellian profile pm{p) ^ 
exp(— p^/ {ImXM))- These distributions will be obtained from 
the numerical data for the oscillators and will be used to check 
the predictions arising from Eqs. (|7]| and ( [T3] l. 

In the limit of large systems, we expect all these mea- 
sures to approach the same equilibrium value; however, for 
small number of degrees of freedom, they are likely to differ 
Il22l . We will use two basic definitions of temperature, and 
the different ways to calculate them, to assess the equilibrium 
properties of our model system. 



DYNAMICS VIA ENSEMBLE AVERAGE (A' = 1) 

In this section, we investigate the approach to equilibrium 
and equipartition of energy for the system described by the 
Hamiltonian ([T]i with only one quaitic system, which means 
N = 1. We examine, particularly, the energy transfer dynam- 
ics of this system for cases in which the coupling between the 
subsystems is weak and the harmonic oscillators are near-to- 
resonance. We test for equilibration by comparing the calcu- 
lated distributions of energy and momentum with appropriate 
equilibrium distributions. We also analyze and compare the 
dynamical behavior of t| with that of T^', concentrating on 
the issue of energy equipartition between the energy stored in 
the oscillators and in the quartic system at large times. 

Our approach makes use of averaging over large ensembles 
of different initial chaotic configurations that have a common 
fixed energy shell where Hqs = Eqs- The initial conditions 
for the two harmonic oscillators are qi = 0, PI = V2mi£i(0), 
q2 = 0, and p2 — VTm^E^iO). The initial data used for the 
phase space variables x,y,p^,py originate from points along 
a single trajectory of the uncoupled chaotic quartic system, 
rather than from random starting points in an energy shell. 

To solve numerically the equations of motion, we used 
the fourth-order Runge-Kutta method (detailed for example 
in 123 il ). The integration time step length was set to ensure 
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Figure 2. Distributions of energy for the three subsystems at time 
tf, computed with initial value Eqs = 5. The ragged black curves in 
each fe-labeled panel are for oscillator 1 ; the gray curves are for oscil- 
lator 2, and the smooth curves are plotted according to Eq. ^16^. The 
bottom-right panel displays results for QS, with k increasing from 
top curve to bottom curve. 
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Figure 3. (Color online) Panel (a): time evolution of (p^) for 
Eqs(O) =0.1 (bottom curves) and Eqs{0) = 5 (top curves). Panel 
(b): time evolution of for Eqs{0) =0.1 (bottom curves) and 
^es(O) = 5 (top curves). Colors and parameters as in Fig. [T] ex- 
cept that, for the calculation of 1^, 4000 initial conditions were used 
instead of 40000. The results for (p^) are nearly identical to those of 



energy conservation to within 1% for each individual trajec- 
tory. The ensemble average value of an observable is calcu- 
lated as the mean of its estimates generated by propagating 
initial conditions in the ensemble. Throughout this section, 
we set mi =m2 — 10 and a = 0.1 iPJl . 

The left panels in Fig.[T]present the short-time energy trans- 
fer dynamics between the two harmonic oscillators. Their 
average energies are plotted versus dimensionless time ~ 
(Oit/ln for the initial energies Eqs = 0.1 (panel (a)) and 
Eqs = 5 (panel (c)), and frequencies ©2 in a small deviation 
from the frequency cOi = 0.0125: CO2 = COi +kA, with fixed 
value A = 0.0000375 for k = {0, 1,6, 11, 16}. The initial en- 
ergies of the oscillators are chosen as Ei = 25 and E2 — 0. 
It is readily seen that the energy curves for k — Q and 1 are 
almost identical for both situations; however the response to 
Eqs = 0.1 is faster than that of Eqs — 5. Another important 
feature is the appreciably less energy transferred from one os- 
cillator to the other at greater k values. 

The right panels of Fig. [Tjin turn capture the long-time be- 
havior of the average energies of the two oscillators under 



Eqs = 0.1 (panel (b)) and Eqs = 5 (panel (d)) conditions. 
Except for the case when the frequencies coincide, the final 
energies for the oscillators are visibly different for each ini- 
tial state of the system. These panels indicate that the equili- 
bration time is very long and might not be reached in practi- 
cal situations. Not shown here is the behavior of the chaotic 
system's average energy, which ascends gradually with time 
and more rapidly with increasing k, approaching a saturation 
value. See panel (a) in Fig.jsjfor a plot of (p^) as a function 
of f*. 

Fig.[2]exhibits the energy distributions of the two oscillators 
for an ensemble of 40000 initial conditions at time tf — which 
corresponds to 3000 periods of the oscillator 1 — for the case 
Eqs{0) — 5. These distributions are generally different (ex- 
cept for A: = 0) and do not have the Boltzmann form (except 
for A: = 1 1 and A: = 16 for oscillator 2). The fittings to these 
distributions are, in majority, of the form 

PfiE)^Ai{ET{tf))-EfE'^, (16) 

where {Ejitf)) is the value of the average total system energy 
at time tf. In most cases, all parameters A, B, and C had non- 
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Figure 4. (Color online) Energy of the oscillators as a function of the scaled time for [(0\,0>i) = (0.3,0.3) (panel (a)), (0.3,0.302) (panel 
(b)), (0.3,0.304) (panel (c)), and (0.3,0.5) (panel (d)). The black and gray solid curves show the numerical results, and black and gray dotted 
curves in panels (a) and (d) the LRT results calculated using Eqs. l |20| l and \2\) with 7~ 0.0044 (panel (a)) and y ~ 0.0035 (panel (d)). 



zero values. One exception occurs for the chaotic environment 
for k equal to zero or one. In such cases, the fitting is reduced 
to A^/E. 

For the case of a single harmonic oscillator, the energy dis- 
tributions should approach a square root line, as suggested in 
fl3J . However, for ^ = 11 and k~\6, the energy distribution 
of oscillator 2 can be well fit by an exponential, which seems 
to be an unexpected transient. Interestingly, the correspond- 
ing momentum distribution closely obey the Maxwellian law 
Pm{p), which also characterizes approximately the momen- 
tum distributions of the quartic system. This behavior per- 
sists for times up to = 6000. In particular, for A; = 11, the 
value found of Tm = 12.28 for QS is reasonably close to three- 
halves of that of Tm = 8.16 for oscillator 2, which in turn is 
near the corresponding value of Tg = 8.77. This satisfies both 
equipartition {EQs{tf)) ~ |t£ and, from energy conservation 
[cf.Eq. ^],TE^^ET{t = 0). 

We show in Fig. dthe time evolution of t| associated with 
(panel (a)) and of (panel (b)) for all k values studied. 
Again, not all curves reach an asymptotic value within the 
displayed time window, confirming that dynamical equilib- 



rium was not yet reached. We compare the asymptotic values 
with those predicted by Eqs. (|7| and ( [T3| l for total energies 

8.57 for panel (a); 



25.1 and 30: f x25.1 



!7.17 and f x 30 ; 



X25.1 



10.04 and fx 30: 



12 for panel (b). Except for the 



resonant cases A: = and k = 1, converges to values close 
to the expected results of 7.17 and 8.57. Also, T^' converges 
better to 12 than it does to 10.04 for the resonant cases. 

Although both and t| display very similar time depen- 
dent behavior in practically all cases studied, they disagree 
with respect to the mean energy at equilibrium. The numeri- 
cal results for the mean oscillators' energies are closer to t| 
than to which indicates that the usual definition of entropy 
[cf Eq. ([12])] can not be applied to this system, possibly be- 
cause of its few degrees of freedom. 



DYNAMICS WITH A SINGLE TRAJECTORY (N = 100) 

In the previous section, we discussed the dynamical behav- 
ior and equilibrium properties of the system described by the 
Hamiltonian ([TJ for N ~ I. The time dependence of any ob- 
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servable, like the energy of the oscillators, typically displays 
large fluctuations when calculated for a single initial condi- 
tion. These fluctuations, which result from the small num- 
ber of degrees of freedom of the environment, are drastically 
reduced when averaged over ensembles of initial conditions. 
Even after such average we cannot state that the environment 
simulates the action of a thermal reservoir, since the energy 
distribution does not always follow a Boltzmann exponential 
law. As more QS's are included in the environment, these sin- 
gle trajectory oscillations decrease. For sufficiently large 
the time behavior obtained for a single trajectory becomes 
similar to that of the ensemble average and independent of A^. 
In this case, we may talk about a potentially effective dynam- 
ics, where single realizations reproduce the average behavior. 

The indirect interaction between the QS's via the harmonic 
oscillators enables the energy to be redistributed among them, 
leading to a Boltzmann type of equilibrium distribution for 
the environment lfT2l . When a single harmonic oscillator is 
in contact with a sufficiently large and "cold" chaotic envi- 
ronment, most of its energy is transferred to the environment. 
The results exhibited in Fig.[T]show that for two resonant os- 
cillators this is not true: more than 50% of the initial energy 
stored in the oscillators remains in the harmonic modes even 
for very long times. In what follows, we explore situations 
where this symmetry is broken for — 100, for which the 
environment is already in the A^-independent regime. All nu- 
merical results in this section are obtained for a single trajec- 
tory and for mj =1112 = I, a = 0.01 ifTSl . and Ei (0) = 10 and 
£2(0)^0. 

Resonant case: Panel (a) of Fig. |4] shows the energy of 
the two oscillators as a function of time for the resonant 
case. Approximately 50% of Ei (0) + E2 (0) remains with the 
oscillators. Notice that the curves displaying Ei (f*) and E2{u) 
are mirror images of each other with respect to the dashed red 
horizontal line at £ = 2.5, which marks the final mean energy 
of the two oscillators. This result is confirmed in panel (a) 
of Fig.[5]where the energy distributions of the oscillators are 
Gaussians centered in £ =2.5. This is by no means a trivial 
result and is in conflict with the equipartition theorem as given 
in Eq. (|7]i, which predicts Ei = £2 ~ 0.072. The condition of 
high symmetry associated with the resonance is responsible 
for this apparent violation of the equipartition theorem as will 
be discussed in the next section. 

Non-Resonant case: Panels (b) to (d) in Fig. |4] show the 
energy of the oscillators for frequencies moving away from 
the resonance. It is still possible to see a flow of energy from 
one oscillator to the other, although this becomes less evident 
as the frequencies become more separated. Also visible is the 
increase in the dissipated energy. Fig. |4] also shows how sen- 
sitive is the dynamical behavior of the system with respect to 
variations in the frequency of the oscillators. In the quasi- 
resonant cases (see panels (b) and (c)), where 0)1/0)2 ~ 1, the 
relaxation time is significantly greater than in the other cases, 
allowing energy exchange between the oscillators for a very 
long time. However, if the ratio 0)1/0)2 deviates more con- 
siderably from unity (panel (d)), the "opacity" of the chaotic 
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Figure 5. Distributions of energy for oscillator 1 (black squares) and 
oscillator 2 (gray circles) for the four cases displayed in Fig.|4] The 
distribution is Gaussian in the resonant case (panel (a), centered in 
E = 2.5), and Boltzmann-like in all off-resonant cases: pb{E) = 
exp{-E/Ei)/Ei with Ei ^ 0.072 and £2 ~ 0.073 (panel(b)); Ei ^ 
0.072 and £2 ~ 0.073 (panel(c)); and Ei ^ 0.075 and E2 ~ 0.069 
(panel(d)). 



medium increases, culminating in almost independent oscilla- 
tors that quickly lose all their energy to the environment. The 
energy distributions corresponding to the cases of Fig. [4] are 
depicted in Fig. [5] and show that the finite environment does 
act as a thermal bath for the two oscillators, both of which 
have the same temperature. Moreover, the equipartition the- 
orem holds true for all three off-resonant cases at long times, 
which means that the relation £1 ss £2 ~ ^Eqs is valid. This 
reinforces the ability of the finite chaotic environment to pro- 
mote dissipation and thermalization |12|. Nevertheless, the 
resonance condition permits an efficient flow of energy be- 
tween the harmonic modes. The aim of next section is to ex- 
plain the special resonant case, and show that in this case the 
equipartition is still valid. 



LINEAR RESPONSE THEORY 

The oscillators obey the dynamics given in the equations 



(17) 



As the dynamics of the variables x„ are chaotic and the interac- 
tion among the QS's is of second order in the coupling, each x„ 
is approximately independent and we may replace X{t) by its 
average {X{t)). Applying then linear response theory (LRT) 
124 J . we get 



iji + Cofq, K-^ ds (pxx {t -s)(qi (s) +q2is)), (18) 
m Jo 
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where the response function (^xz is given in Ref. lfT2l by 



4 ds 



dtds 



(19) 

where is a parameter that depends on the average energy of 
the environment. After computing the integral we have 



(20) 
(21) 



with 7 = . These equations cannot be diagonalized unless 
Oi = 0)2. In this case, we can define the new variables 



qi ±q2 



(22) 



associated with center-of-mass and relative coordinates of the 
oscillators, and rewrite ( |20] i and ( |2T] i as 



64 



2724 

Q- 



0, 



co"e_ =0. 



(23) 
(24) 



The 2- is then completely decoupled from the environment, 
and will be called the conservative mode, whose initial en- 
ergy is £'i(0)/2 for the present initial conditions. Analo- 
gously, 2+ will be termed dissipative mode. LRT predicts 
7^0); therefore Q+ dissipates energy according to £+(?) ~ 
£'+(0)exp(-27f), where £+(0) is also £'i(0)/2. Note that 
Eq. p4l ) is exact, as can be seen by replacing ( |22] l into ( fTTj i. 

The harmonic oscillators energies obtained from Eqs. ( pOj i 
and ( |2T] l are compared in Fig.[4]with numerical simulations of 
the Hamiltonian ([TJ. As seen from the figure, the qualitative 
agreement is excellent for both resonant (panel (a)) and non- 
resonant (panel (d)) cases. 
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Figure 6. (Color online) Distribution of energy for the dissipative 
2+ -mode for A' = 100. The fit is pb{E) = exp (-£/£+)/£+ with 
E+ = 0.037. 

In order to correctly apply the equipartition theorem in 
the resonant case we have to consider only the energy in 



Q+. Taking the expression for the total initial energy, Ef — 
£1 /2 + NEqs, we find that the expected value of the energy 
of the dissipative mode in equilibrium is E+ w 0.039. This 
is confirmed in Fig. [6] which shows the energy distribution 
of Q+. Therefore, the violation of equipartition was only ap- 
parent, caused by the emergence of a conservative mode that 
prevented the complete dissipation of the energy in the har- 
monic oscillators. 



CONCLUSIONS 

We have studied the dynamical behavior of two harmonic 
oscillators independently coupled to a chaotic environment 
with 2N degrees of freedom. We focused our analysis on two 
points: the interactions between the oscillators mediated by 
the chaotic environment and the equilibrium properties of the 
system. We found that the oscillators can exchange energy 
through the environment when in almost perfect resonance. 
Deviations from this condition quickly changes this behavior 
and makes the oscillators less sensitive to the presence of each 
other When in perfect resonance, the oscillators are able to 
keep part of their initial energy in an apparent violation of the 
equipartion theorem. This, however, turns out to be exactly 
the fraction of the energy stored in the oscillators relative mo- 
tion, which is not coupled to the environment. This holds both 
for = 1 via ensemble averages or for = 100 for a single 
realization of the dynamics and is also true within the approx- 
imation of the linear response theory, which works well for 
the resonant and non-resonant cases. 

The equilibrium properties of the system, on the other 
hand, depend critically on the number of degrees of free- 
dom of the environment. In order to quantify the equilibrium 
we considered two measures of temperature: T£, obtained 
from the equipartition theorem, and T5, obtained via the usual 
definition of entropy. These temperatures were also com- 
pared with the energy distribution of each subsystem, when- 
ever these converged to a Boltzmann-like exponential. For 

= 100 the oscillator's energy distributions indeed always 
converged to the exponential decay exp (— £ / Tg) with Tg w Te 
for all subsystems, except in the resonant cases. This charac- 
terizes the thermal equilibrium and corroborates earlier results 
that a not too small chaotic environment does play the role of 
a thermal bath 1, 12J . This conclusion also holds in the reso- 
nant case if the energy of the conservative mode is properly 
subtracted. 

The case A^ = 1, on the other hand, showed a very rich be- 
havior Equilibration takes very long times for the present 
choice of parameters and the energy distribution of the oscilla- 
tors display, in some cases, exponential curves. Curiously, for 
the particular values of coi and 0)2 used, the exponent agreed 
reasonably with energy equipartition, but not with Tg- More- 
over, the distribution of the momentum did often displayed 
the expected Gaussian distribution, even if the corresponding 
energy distribution was not exponential. These unexpected re- 
sults show that the N = I environment does not simulate the 
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action of a thermal bath even if averaged over an ensemble of 
trajectories, and display a more complicated type of approach 
to equilibrium that is worth a deeper investigation. 

We have also examined the ergodicity of the system. Be- 
cause X is small, it could be expected that the total system 
would not be ergodic. This, however, is not the case. For 
A = 0, the environment is ergodic for N = 1 but it is not er- 
godic for — 100, since the quartic systems are totally inde- 
pendent from each other As the coupling is turned on, the to- 
tal system becomes ergodic for = 100 (and is able to show 
the correct equipartition of energy), whereas, for = 1, it 
does not. This was verified by numerical results not shown in 
which time and ensemble averages were compared for various 
conditions. As expected, we found that temperature could be 
defined whenever ergodicity was satisfied, although ergodic- 
ity itself depends non-trivially on the parameters of the system 
and on the number of degrees of freedom. 

Finally, we note that the emergence of the conservative 
mode Q depends on the symmetry of the coupling as given 
in Eq. Q and on the resonance condition coi = (02- The con- 
servative mode is absent if Oi 7^ 0)2 or if each oscillator is 
coupled to a different mode of the quartic system, such as in 
Hi = Y^nili^n +<j2yn)- In LRT this leads to decoupled equa- 
tions for 2+ and each equation identical to Eq. ( |23| ). 
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